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Abstract 

The effect of an external pressure disturbance, being displaced with a constant speed along the free surface of a 
viscous thin film, is studied theoretically in the lubrication approximation in one- and two-dimensional geometries. 
In the comoving frame, the imposed pressure field creates a stationary deformation of the interface - a wake - that 
spatially vanishes in the far region. The shape of the wake and the way it vanishes depend on both the speed 
and size of the external source and the properties of the film. The wave resistance, namely the force that has to 
be externally furnished in order to maintain the wake, is analyzed in detail. For finite-size pressure disturbances, 
it increases with the speed, up to a certain transition value above which a monotonic decrease occurs. The role 
of the horizontal extent of the pressure field is studied as well, revealing that for a smaller disturbance the latter 
transition occurs at a higher speed. Eventually, for a Dirac pressure source, the wave resistance either saturates for 
a ID geometry, or diverges for a 2D geometry. 


1 Introduction 

A disturbance moving atop a fluid or a soft medium deforms the interface shape of the latter, thereby producing a 
wake. This is a common phenomenon, which takes place in various natural, scientific and industrial settings [T] across 
several orders of magnitude, ranging from waves generated at the surface of water by animals and watercrafts [2], to 
patterns observed during the deposition of thin polymer coatings [3]. For instance, at the geophysical scale, the study of 
gravity-driven flows around human-made obstacles and landscape singularities is essential to understand and minimize 
the damages caused by natural disasters, such as snow avalanches [4], landslides [5] and lava or mud flows laiiiiH]. 
In nanophysics, gaining insight on the intrusive effect of local probes, such as atomic force microscopy [9l[T0l[IT], is 
crucial in order to improve further the measurement accuracy. 

Similarly, the wake angle and the wave pattern generated behind a disturbance sailing at the surface of an inviscid 
liquid remain current topics of interest na Hal nn, because they are crucial for naval industry. In fact, to be 
maintained, the wake continually consumes energy, the latter being radiated away from the disturbance. This loss can 
be formulated into a force, called the wave resistance m, which needs to be furnished by the controller in order to 
ensure constant speed. The wave resistance associated to the capillary-gravity wakes has thus been fully characterized 
for an incompressible inviscid flow HU nail], including the cases of Dirac and finite-size pressure distributions. The 
perturbative role of viscosity on the inertial wave resistance m has also been considered, as well as the effect of 
thickness for a viscous liquid film BO]. 

In the context of viscous thin films Eauana, a moving external disturbance and the associated wake at the 
free surface may directly be used as a new kind of fine rheological probe |24l [25] . This idea was also proposed as a 
possible method to measure slip at the solid-liquid boundary [26] . Interestingly, despite the broad applicability of the 
lubrication theory and its importance for industrial processes, the study of the wake and the wave resistance in such 
a context is a topic which has been hastily explored and remains an open question. For instance, the competition 
between visco-capillary relaxation and perturbation speed, which might control the extent and lifespan of the wake 
generated by the moving perturbation, has not been studied. Imagine that the footprint left behind by the moving 
perturbation is arrested by performing a local quenching, before the relaxation of the thin-film occurs Emuj; this 
may allow to imprint a narrow channel in the film, thus revealing the foundations of a new nano-patterning technique. 
Additionally, one may be interested in estimating and controlling the amount of energy that this technique requires, 
for which the wave resistance of a moving disturbance atop a viscous thin film must be well characterized in detail. 

In the present study, the surface pattern generated by a moving external pressure disturbance along a viscous thin 
film is described theoretically, within the lubrication approximation, for one- and two-dimensional systems. The wave 
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2 GENERAL CASE 


resistance is analyzed in both cases, leading to analytical expressions of this force in terms of the disturbance speed 
and size, and the physical properties of the liquid film. Low- and high-speed regimes are described through asymptotic 
expressions, which unveil the nature of the wave resistance in the lubrication context. 


2 General case 

2.1 Thin-film equation for a moving pressure disturbance 

We consider a viscous liquid film of thickness ho deposited over a flat horizontal substrate. The equilibrium free 
surface of the film is flat as well, and perpendicular to the vertical direction 2 ;. An external pressure field -0, moving 
along the horizontal direction x with constant speed v > 0, is applied over the liquid film. As a result, the incom¬ 
pressible Newtonian fluid flow yields a stationary surface pattern h (see Fig.[^ in the frame of reference of the moving 
disturbance. 

In the lubrication approximation, neglecting shear stress at the liquid-gas interface and assuming no slip at the 
substrate, the incompressible Stokes equation and conservation of volume lead to the following Reynolds equation [ 21 ]: 


dh _ 

m=^- 


{ho + hy 

3/i 


Vpig 


( 1 ) 


where fi is the dynamic viscosity of the liquid and pig is the pressure jump at the liquid-gas interface. Under the 
lubrication hypothesis for the no-slip case, one can also neglect the normal viscous stress at the liquid-gas interface, 
and the pressure jump is thus given by the Young-Laplace equation. For small slopes, the latter reads: 


Pig = -yAh - 1 - pgh -h y 


( 2 ) 


where A denotes the Laplacian operator in cartesian coordinates, g is the acceleration of gravity, and 7 and p are the 
liquid-gas surface tension and density difference, respectively. Equations 0 and 0 are coupled, and should be solved 
together in order to describe the dynamic response of the viscous thin film due to the displacement of the external 
pressure field. 

We introduce the dimensionless variables: 
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where the capillary length k~^ = \/'y/{pg) has been chosen as the characteristic length scale in the a;-direction, and i is 
a reference length scale in the p-direction, which depends on the dimension of the system: in ID one sets £ = Ly, a unit 
length, and in 2D one sets i = In turn, r = Spy/ [{pg^ho] is the visco-capillary time scale, and Sext = f f dxdyy 
is the imposed load, where in ID the integral over y is restricted to the horizontal extent Ly. In the limit of small 
surface deformations H 1 , the thin-film equation (TFE), obtained from the combination of eqs. 0 and Q, is given 
by: 

— = -V ■ [V (Aff) -VH- v^] . (4) 

Let us introduce the capillary number Ca = pv/^ and the Bond number Bo = pga‘^1^, where a denotes the 
characteristic horizontal size of the external pressure field. We now place ourselves in the frame of reference of the 
moving disturbance, through the new variable U = X — UT, where the reduced speed V reads: 


V = VTK . 


(5) 


By assuming stationarity of the surface profile in this comoving frame, one has H (X, Y^T) = ( (U^Y) and thus eq. 0 
reduces to the partial differential equation: 
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In the following, we refer to eq. 0 as the 2D-TEE. 

Similarly, for a ID geometry, in which the system is invariant with respect to the E-direction, eq. 0 is reduced 
to the ordinary differential equation: 
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2.2 Wave resistance 


3 ONE-DIMENSIONAL CASE 



Figure 1 Schematic diagram of the surface profile h of a viscous thin film reacting to an external pressure disturbance L 
moving along time t at constant speed v. (a) ID and (b) 2D geometries. 


After integrating once with respect to U, and considering the boundary conditions C = 0? d(/dU = 0, d^(/dU^ = 0, 
and d^/d[/ = 0 at \U\ oc - since the liquid-gas interface should be undisturbed far from the moving disturbance, 
and since the flow and pressure gradients vanish in the far field - one gets: 

df/3 dU 

In the following, we refer to eq. as the ID-TFE. 


2.2 Wave resistance 


The wave resistance r is the force that has to be externally furnished by the operator - i.e. the pressure field - in 
order to maintain the wake. According to Havelock’s formula m, it reads: 


dh 




(9) 


where in the ID case the integral over y is restricted to the chosen horizontal extent Ly. Defining the dimensionless 
wave resistance R through: 


R = 


and recalling eq. (§, one gets in the comoving frame: 
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( 10 ) 

( 11 ) 


In the next two sections, we solve the ID-TFE and 2D-TFE using Fourier analysis, and study the shape of their 
solutions and the wave resistance for different values of the two relevant dimensionless parameters: the reduced speed 
V and the Bond number Bo. One should keep in mind that ^/^o = an is the rescaled lateral extent of the external 
pressure field. Also, note that, by linearity of the response, all the results correspond to an arbitrary third dimensionless 
number S'ext, the dimensionless load, as it can be absorbed in C,. In the following, S'ext is chosen to be positive, although 
the physical foundations and the wave resistance are valid for any sign of S'ext- 


3 One-dimensional case 

3.1 Surface profile 

Considering the one-dimensional Fourier transforms ^ (Q) and (Q) (see App. of the surface profile and the 
pressure field, the ID-TFE given in eq. (§ becomes: 

— {iQ^ + iQ + F) C = iQ^Sext^ ? (12) 
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3.1 Surface profile 


3 ONE-DIMENSIONAL CASE 
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Figure 2 Normalized surface profile solution of eq. Q as given by eq. (13), using Lorentzian (see eq. (14)) and Dirac 

{Bo 0) ID pressure distributions, with dimensionless load iSext > 0, for different values of the reduced speed V and the 
Bond number Bo^ as indicated. 
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3.1 Surface profile 


3 ONE-DIMENSIONAL CASE 



V 


Figure 3 Real (plain) and imaginary (dashed) parts of the three dimensionless wavelengths Xj of the asymptotic solution, as 
given by eq. (17) for the ID case, as a function of the speed ratio V = V/Vc = 3\/3y/2, where V is the reduced speed. Note 
that for V < 1, (Aj) = 0 for all the values of j, whereas for V > 1, Re {X 2 ) and Re (A 3 ) overlap. 


where i = is the imaginary unit, and Q is the angular wavenumber in the [/-direction, 

profile C is described by: 


C{U) 


^ext f (Q) exp {illQ) 

^ J iV-Q{Q^^l) ^ ’ 


As a consequence, the 


(13) 


where we recover the statement above that the dimensionless load Rext gives only a constant prefactor in the linear 
response. 

In order to account for the finite-size effects, we introduce a Lorentzian pressure distribution: 


^{U) = 




Tr{m +Bo) ’ 

(Q) = exp (-\/^|Q|) . 


(14) 


Note that f dU^ = 1, by construction, and that in the limit Bo 0 this pressure distribution becomes ^ (U) = S ([/), 
or ^ (Q) = 1, where 6 denotes the Dirac distribution. 

The associated normalized surface profile is shown in Fig. for different values of the reduced speed V and the 
Bond number Bo. At low speed, the deformation profile is found to be nearly symmetric. In that case, the depression 
is centred at the maximum of the pressure field and it exponentially relaxes towards a fiat horizontal surface, far away 
from the perturbation. As we increase the speed V, the symmetry of the surface shape breaks and the depression is 

shifted towards the rear (left-hand side of Fig. of the driving pressure field, while a bump at the front (right-hand 

side of Fig.|^ grows with the disturbance speed. An exponential relaxation is still observed in the rear, with a shorter 
characteristic length scale as the speed increases, whilst either a similar relaxation (for intermediate speeds) or an 
exponentially decaying oscillation (for high speeds) is discerned at the front. 

These features can be understood by looking at the asymptotic solution of the ID-TFE. At a sufficiently far position 
{\U\ :$> \/Bo) from the pressure field, its overall effect can thus be taken as an interior boundary condition in order to 
analyse the spatial spectrum of the solution in the far field. In that approximation, eq. § has the three asymptotic 
basis solutions: 

Ci ~ exp (^E^ ^ ( 15 ) 

for j = 1,2,3, where the wavelengths Xj are given by the characteristic equation: 

“ 1) = ’ ( 16 ) 
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3 ONE-DIMENSIONAL CASE 


and thus: 
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Here, V is the speed ratio V = V/Vc defined from the critical speed Vc = 2/(3\/3). The dependence of \j on V is 
depicted in Fig. Since Q comes from an expression which may have different values (the roots of — 2VO^ + 1) 
we note that the correct value to be employed in eq. ( [Tt] ) satisfies the relation [H = 1 — [^ (^)]^ for V < 1 , and 
= 0 for V > 1 . 

According to the signs of Xj in Fig. and in order to avoid any divergence of the profile at U = — oo, it follows 
that for U <0 and |t/| ^ \/Bo the shape of the surface is necessarily given by: 


( = Ni exp 



(18) 


where Ni is a real constant and Ai takes a real value 5)?(Ai) > 0. This indicates that the surface always presents a 
pure exponential relaxation towards a fiat horizontal profile at the rear of the perturbation, in the far field, no matter 
the speed. The typical distance Ai over which this relaxation occurs becomes shorter as V increases, as observed in 

Fig.[l 

On the other hand, for 1/ > 0 and |t/| ^ \/^, the behaviour of the asymptotic solution depends strongly on V. 
When V < 1, according to the signs of Xj in Fig. and in order to avoid any divergence of the profile at U = +oc, 
the surface profile is described in the far field by the superposition: 


c 


N 2 exp 



+ Ns exp 



(19) 


where N 2 and Ns are real constants, and A 2 and A 3 are both real and negative with IA 2 I < jAsj. As a consequence, 
an over-damping of the surface perturbation is observed in the far field. The V = 1 case corresponds to the critical 
damping situation for which A 2 = A 3 = — 1. Finally, for V > 1, A 2 and A 3 are complex conjugates, and so are the 
corresponding amplitudes N 2 and Ns. Thus, one can write: 


C = N exp 


^{X2)U 
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( 20 ) 


where N and 0 are real constants. This indicates the presence of spatial oscillations at the front of the perturbation in 
the far field, before the film surface attains its complete relaxation. In brief, the critical speed W denotes the reduced 
speed V at which a transition is observed from overdamped relaxation to damped oscillations, at the front of the 
perturbation. 

In addition to the far-held trends described above, that depend solely on the reduced speed V and the associated 
dimensionless characteristic lengths Aj, one can discern in Fig. [^different near-held behaviours of the surface prohles 
that involve a second dimensionless length ^/Bo associated with the horizontal extent of the pressure held. As already 
explained above, for a pressure distribution with sufficiently small size (\/^ ^ 1 ), the healing length of the prohle is 
almost independent of Ho, and approximately given by Aj that decreases as V increases. In contrast, for a large size 
^ 1), the lateral extent of the prohle is comparable with the typical size \/Ho of the pressure disturbance, as 
can be observed at low speed V in Fig. In this ^/Bo ^ 1 case, the healing length decreases with increasing speed 
only for low velocities. Above the critical speed W typically, the healing length saturates, leading to surface prohles 
that nearly follow the same trend, independently of V. 


3.2 Wave resistance 


Invoking eq. (11) and the fact that in ID one has i = Ly^ one gets the following ID expression of the dimensionless 

R = 


wave resistance: 


5ext J {U) ^dU . 


( 21 ) 


For a viscous thin him, it can be shown that the power exerted by the wave resistance equals the viscous dissipated 
power in the bulk (see App. eq. (58)). This power balance highlights the fact that the wave resistance is the force 
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3.2 Wave resistance 


3 ONE-DIMENSIONAL CASE 



log(l/) 


Figure 4 Dimensionless wave resistance R (see eq. ( |22[ )) normalized by the square of the dimensionless load S^ext, as a function 
of the reduced speed V, for a ID Dirac pressure distribution, and for a ID Lorentzian pressure distribution (see eq. (14)) with 
different values of the Bond number Bo as indicated. The vertical dashed lines indicate the crossover speeds Vs, estimated 
from eq. (30) for different Bo. The horizontal dotted line corresponds to the plateau value of 1/6 in the ID Dirac case. 



log(B,) 


Figure 5 

Bo. Both 


Shape function / (see eq. (23)) of the ID pressure field in the low-speed regime, as a function of the Bond number 
the Lorentzian (plain, see eq. (24)) and the Dirac (dashed, see text) cases are represented. 
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4.1 Surface proRle 


4 TWO-DIMENSIONAL CASE 


externally furnished by the operator in order to maintain the stationary wake. Remarkably, due to the nature of the 
TFE, the wave resistance is also related to the spatial energy spectral density (see App. eq. ([4^). 

Combining eqs. (13) and (21) leads to: 


R = 


SLtV f Q2$(q)§(_q) 


/ 


2^ J + 


dQ . 


( 22 ) 


Figure [^displays the normalized wave resistance for a Lorentzian pressure distribution, as given by eqs. (14) and (22), 
as a function of the reduced speed F, and for different values of the Bond number Bo. The wave resistance for a Dirac 
pressure distribution is also shown. 


For the asymptotic case where F ^ 0, a series expansion of the integrand in eq. (22) leads to: 


R = 


[f{Bo) + 0{V^)] , 


(23) 


where / {Bo) is a function that depends only on the shape of the pressure field. This points out the linear dependence 
of the wave resistance on F, in the small speed regime, as observed in Fig. For a Dirac pressure distribution one 
gets fs {Bo) = 7r/4, whereas for a Lorentzian (see eq. @) , one has: 


pc 

/l {Bo) = 8Bo3/2 / 
Jo 


exp {—w) 
{w^ + 4Bo)" 


:dw 


(24) 


Both behaviours are depicted in Fig. In addition, /l ^ 7r/4 when Bo ^ 0, as for the Dirac pressure distribution, 
whilst a series expansion at Bo oc yields: 


/l {Bo) = 


2v^ 


[l + 0(5o-i)] , 


(25) 


observed in Fig. 

Considering the Dirac delta distribution, for any F, one finds (see App. [A|for details): 


Rs 


-iV30^/^Sl,V 

2 (6>3 - 26>2 + 29-1) ’ 


with: 


6> = Y1 - 2V2 + 2vyv2 - 1 , 

where V was defined after eq. (0. For the asymptotic case where F ^ oc, a series expansion leads to: 


Rs = ^ 

D 


1 + 0 




(26) 

(27) 

(28) 


This result indicates the singular feature that for a ID Dirac pressure distribution the wave resistance saturates to a 
constant value at high speed, as observed in Fig. 

For a finite-size Lorentzian pressure distribution, we combine eqs. (14) and ( [^ (see App. for details on the 
overall expression) and obtain in the high-speed regime F ^ oo: 




(29) 


Therefore, a 1/F dependency of the wave resistance is asymptotically expected for finite-size pressure fields at large 
speed, as indeed observed in Fig. 

For a Dirac pressure distribution, the reduced speed Vs at which the wave resistance saturates can be defined 
by equating the low- and high-speed asymptotic expressions of Rs. This provides F = 2/3, consistent with the 
observation of Fig. Similarly, for a finite-size Lorentzian pressure distribution, the reduced speed F at which the 
crossover between the low- and high-speed regimes of Rl occurs is defined by: 


Fs 


2Bo3/4y7(^ 


(30) 


as observed in Fig. In the case of small Bond number Bo, I 4 scales as Bo 3/^; whereas for large values of Bo, I 4 
scales as Bo~^^^ . 
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4.1 Surface proRle 


4 TWO-DIMENSIONAL CASE 




Figure 6 Normalized thin-film surface profile in presence of a moving Lorentzian pressure disturbance, with dimensionless 
load iSext > 0, Bond number Bo = 10^, and reduced speed V = 10, computed from eqs. (32) and (33). Perspective views of the 
front (top) and the back (bottom), in arbitrary units and colour code. 
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Figure 7 Normalized thin-film surface profile (top view) computed from eqs. (32) and (33), for a dimensionless load Se^t > 0, 
and different values of the Bond number Bo and the reduced speed V. In each subfigure, the pressure field travels from left to 
right, along the longitudinal coordinate U. 
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4 TWO-DIMENSIONAL CASE 


4 Two-dimensional case 


4.1 Surface profile 

Considering the two-dimensional Fourier transforms ( (Q, K) and ^ (Q, K) of the surface profile and the pressure field, 
with the definitions given in App. the 2D-TFE of eq. §) becomes: 


iQV - (Q2 +K^)^- (Q2 + (g2 + 7^2) 5 ^ 


(31) 


where Q and K are the spatial angular frequencies in the U and Y -directions, respectively. Consequently, the surface 
^ {U, Y) of the thin film reads: 


< ii iQV-(Q^+K^)(lYQ^YK^) ' 


(32) 


In analogy with the previous ID case, we introduce the following 2D axisymmetric finite-size Lorentzian pressure 
distribution: 

y/Bo 


^{U,Y) = - , 

27r (f/2 -h y2 + 

^{Q,K) = exp [- ^/Bo{Q^ + K^) 


(33) 


We recall that, since a is the characteristic horizontal size of the pressure field, \/Bo = an is the dimensionless size. 
When, Bo 0, the Dirac pressure distribution ^ {U^Y) = 6{U)6{Y), or its equivalent ^ {Q^K) = 1 in the frequency 
domain, is recovered. 

Surface patterns generated by a Lorentzian pressure disturbance, for several combinations of the parameters Bo 
and E, are illustrated in Figs. and A common feature of these profiles is the depression observed. In Fig. due 
to the particular parameters employed (see caption) to generate the surface profile, a non-symmetric shape is shown. 
A rim surrounds the front of the depression and becomes a wake, behind, that follows the straight translation path of 
the pressure distribution along U. 

In Fig.[^ for low speeds {V < 1), we find deformation profiles which are mostly symmetric, with a main depression 
aligned with the maximum of the pressure field and followed by a monotonic relaxation towards a flat horizontal 
surface. In this situation, the liquid-film surface relaxation occurs faster than the displacement of the pressure field: 
in real units, the displacement speed v is smaller than the characteristic speed pghQK/3jii at which the film relaxes and 
thus propagates information at its surface. As the pressure acts on the surface by creating a moving depression, the 
film has enough time to move the dislodged liquid volume isotropically. 

Like for the ID case, an increase of the speed V induces a symmetry breaking phenomenon, which above V 1 
consists of a level rising at the front (right-hand side of the corresponding sub figures of Fig. and a wake at the rear 
(left-hand side). In real units, the relaxation of the liquid film now occurs with a speed which is smaller 

than the displacement speed v of the pressure field. In other words, the liquid volume that is ejected by the applied 
pressure is not distributed isotropically because the information at the film surface does not propagate fast enough. 
This retardation effect is similar to the Mach and Cerenkov physics [T], where the speed of the object becomes larger 
than the wave-propagation speed in the considered medium. 

For relatively moderate speeds, E ^ 1, the protuberance at the front covers a crescent-shaped region that grows 
with the size of the pressure field y/Bo. The maximum surface level at the front of the pressure field is located at 
T = 0, whereas the maximum of the wake, at each position U < 0, is located at T ^ In contrast, the local 

minimum of the wake is always placed at T = 0. In addition, both the local maximum and minimum of the wake 
present a rapid spatial exponential relaxation. As the speed grows in the range E G [l, 10^], the height of the frontal 
rim reduces, while the extent of the wake in the rear reaches a larger distance. The local maximum at the front is still 
located at T = 0, while the maximum height of the wake is now placed at Y ^ \U\^^ where the exponent a evolves 
from 1/2 at E to 1/4 at E ^ 10^. In fact, the 1/4 exponent in the high-speed regime indicates the same dynamics 
as the capillary leveling of a trench in a viscous thin film with no slip at the substrate EH EH]. In this regime, the film 
seems to overlook the presence of the fast disturbance, perceiving the footprint left by the latter as a ID trench which 
the film tends to refill by capillary leveling. Therefore, when we compare the transverse profiles at different positions 
U < 0, we in fact discern different stages of the ID leveling of a viscous film. 


4.2 Wave resistance 


Invoking eq. and the fact that in 2D one has £ = k, one gets the following 2D expression of the dimensionless 
wave resistance: 


R = S,^tjj^ (V Y) ^dUdY . 


(34) 
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4.2 Wave resistance 


4 TWO-DIMENSIONAL CASE 



log(V) 


Figure 8 Dimensionless wave resistance R normalized by the square of the dimensionless load Sext, as a function of reduced 
speed y, for a 2D Dirac pressure distribution (see eq. (62)), and for a 2D Lorentzian pressure distribution (see eqs. (33) and 
(36)) with different values of the Bond number Bo as indicated. The vertical dashed lines indicate the crossover speeds Vs, 
estimated from eq. (42) for different Bo. 



log(Bo) 


Figure 9 Shape function F (see eq. (37)) of the 2D pressure field in the low-speed regime, as a function of the Bond number 
Bo. Both the Lorentzian (plain, see eq. (38)) and the Dirac (dashed, see text) cases are represented. 
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5 CONCLUSION 


Therefore, in combination with eq. (32), one obtains: 

f f (Q" + K^) $ (Q, K) $ (-Q, -K) 


R = 


"^ext 

47r2 


// 


Q2y2 + (Q2 + J^2f (1 + Q2 + ^^2)2 




(35) 


Considering the polar coordinates (^, (p), an axisymmetric pressure field ^ (p), as it is the case for the Dirac and 
Lorentzian distributions presented in this study, and performing an integration over ip (see App. H yields: 


02 

D ^ext 


f 


1 - 


Q + l) 


^V^ + Q^Q-^ + l)\ 




1 2 


g^dg . 


(36) 


For a Lorentzian distribution, a combination of eqs. (33) an d ([M| ) leads to the corresponding wave resistance, which 
behaviour in terms of the reduced speed V is shown in Fig. for different values of Bo. The wave resistance for a 
Dirac distribution is also represented. 

In the asymptotic case where y ^ 0, a series expansion of the integrand in eq. ( [3^ (see AppJ^for details) leads 
to: 

R=^^[FiBo) + 0{V^)] , (37) 

where F (Bo) is a function that depends only on the shape of the pressure field, pointing out the linear dependency 
of the wave resistance on V, in the low reduced-speed regime. For a Dirac pressure field Fs {Bo) = 1, whereas for a 
Lorentzian, it is described by: 

exp {—w) w 


no 

Fi^ {Bo) = SBo 

Jo 


-dw . 


(38) 


/o (w2 + 4Bo)2 

Both trends are illustrated in Fig. For the Lorentzian Fl {Bo) = 1 at Bo 1, the same behavior as for the Dirac 
pressure distribution, while a series expansion at Bo oc gives: 


1 


F^(Bo) = —J1 + 0{Bo-^)] 


Coming back to the Dirac pressure distribution, a series expansion of eq. (36) at F ^ oc yields: 

F|Mil52^^yi/3[i + 0(y-i)] , 


Rr = 


which indicates that the wave resistance of a 2D Dirac pressure diverges slowly with increasing V. 

For a Lorentzian pressure distribution, a series expansion in the high-speed regime where V 
reduce the combination of eqs. (36) and (33) into: 


(39) 


(40) 


oo allows us to 


Rl = 


35elt [l + 0(y-2)] , 


167rVBo‘^ 


(41) 

OO, the wave resistance 


This demonstrates that, for finite-size pressure fields in the high-speed regime where V 
behaves in the same way as in the ID case, i.e. decreasing as 1/V. 

For a Dirac distribution, the speed dependence of the wave resistance changes from R r\j V to R^ FV3 at the 
reduced speed I 4 = 5/3, which is obtained by equating the low- and high-speed asymptotic expressions of Rg. Similarly, 
for a finite-size pressure distribution, the reduced speed V^, at which the crossover between the two asymptotic regimes 
takes place, reads: 


"" y 2Bo‘^F (Bo) ’ 

emerging when we match the low- and high-speed regimes of Fl- At small Bond number Fo, the crossover speed is 
approximately given by ^ ^Fo^2/3^ , whilst for large values of Fo, it happens around Vs r\j (Bo/3)-'/l 


5 Conclusion 

In the present work, the effects of a pressure disturbance moving near the free surface of a viscous thin film have been 
studied in detail. Depending on the disturbance-speed regime, the existence of nearly symmetric and non-symmetric 
profile shapes have been revealed by means of numerical methods. The nature of the wakes and bumps observed at the 
surface, as well as the lifespan of the trail and its relaxation towards a flat surface in the far field, have been explained 
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through asymptotic solutions, unveiling their reliance on several dimensionless parameters: the Bond number Bo^ the 
dimensionless load 5'ext and the dimensionless speed of the disturbance V, which are related to the geometrical and 
physical properties of the liquid film: the density p, the air-liquid surface tension 7 , the dynamic viscosity p, the film 
thickness Hq and the pressure distribution through its total load 5ext and width a. 

In addition, the wave resistance, which one must overcome in order to displace the disturbance at constant speed, 
has been analyzed. Even though the shape of the free surface is calculated numerically, the wave resistance can 
be derived in different asymptotic speed regimes, by means of analytical expressions for Dirac delta and Lorentzian 
pressure distributions. Moreover, a general expression of the wave resistance for a Dirac delta distribution in the ID 
case has been presented. 

We believe that the measurement of the wave resistance can be employed to deduce one of the physical properties 
of a viscous thin film, such as the above p, 7 , p or ho, when the remaining parameters are known. The working 
principle of a new fine rheological probe may be based on this acquired insight. Additionally, in combination with 
local quenching, controlling the lifespan of the footprint generated by the moving perturbation might serve as the 
technological foundation of a new nano-patterning technique, which may allow for imprinting nano-channels over 
polymeric materials. This idea may provide an alternative to the current nano-lithography techniques employed to 
fabricate nano-devices and templates for electronics. 
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A Details on the one-dimensional case 

Fourier transform 


r° 

f{Q)= / 

J —C 


/ (U) exp i-iUQ) dU , 


f(U)= — 


1 r 

2^ J-, 


f (Q) exp {iUQ) dQ . 


(43) 


Wave resistance 

Multiplying eq. § by the surface displacement, and integrating over the whole U domain gives: 


(44) 


Some of the terms (the first two from the left-hand side and the one in the right-hand side) of this equation can be simplified to: 


/ 


d^C 


dU^ 


I 

/ 


C — dU = 
dU 

C^du = 
dU 




C'l' 


2 \dUJ 


= 0 , 


= 0 , 


- oo 

f 'if — dU = - 

f 'if — dU 

- —oo 

1 dU 

1 dU 


(45) 


if one takes into account that for \U\ —)■ oo, the surface has completely recovered its flat shape, and thus that ^ and all its derivatives are 
equal to zero in the far field. Therefore, one finds: 


Sext J = V j C^dU , 


(46) 


where the left-hand side term is the definition of the wave resistance given in eq. and the integral in the right-hand side is the spatial 
energy spectral density. _ 

Injecting the Dirac delta pressure distribution into eq. \22\ and performing the integration, yields the following expression of the wave 
resistance: 


Rs = 


2 (^1 + q 2 ) (<7l + qz) {q 2 + qs) 

where i = \/—l is the imaginary unit, and where we introduced the auxiliary quantities qi, q 2 and q^, which are given by (see Fig. |l0| >: 


(47) 


' = (l - 2V2 + 2V\/V2 - ij 


1/3 


(48) 


and where the (jj and V have been introduced in eq. After som e algebra, one is able to regain eq. ( |26| >. 

From the substitution of a Lorentzian distribution into eq. \22\ and after integration, the following expression of the wave resistance 
is obtained: 


i^L = 


-^ext _ Qj 

27r j=-s ^ ^ ^ 


f 


— exp {—w) 
w — 2Qj \/Bo 


dw , 


where Qj is the jth root of (Q^ + given by: 


Vs \o-. 


6 


1/2 


(49) 


(50) 


for j = 1,2,3, whereas for j = 4,5,6 we have Qj = —Qjs. (T i an d 6 have been introduced in eqs. 0 and ( |48| >, respectively. The 


asymptotic behaviours for V —)■ oo can also be obtained from eq. 


. 0, 


after making 0 —)■ 1 and some algebra. 
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V 


Figure 10 

speed ratio 
overlap. 


Real (plain) and imaginary (dashed) parts of the three dimensionless auxiliary quantities qj as a function of the 
V, as given by eq. (48). Note that for V < 1, 3^ (g^) = 0 for all the values of j, whereas for V > 1, ^ (^i) and ^ (^ 3 ) 


Viscous dissipation 

The viscously-dissipated power indicates the rate at which the external work is irreversibly converted into internal energy by viscous forces. 
For a thin film, the viscously-dissipated power per unit length is expressed as: 




IH: 


ho 


dvx 

dz 


dz 


dx , 


(51) 


where Vx is the horizontal component of the velocity within the film. Under the assumptions of the lubrication theory, the previous 
expression can be re-written as: 

r u.dApY 

O-h)- 


IJiL 


dx 


dz > dx . 


(52) 


Integration over the local film thickness and, afterwards, integration by parts yields: 

2 

dx , 




|(fto + ft)^(^) 




= I 


(ho + h)^ f dAp\ 


Sfi \ dx J 


dx 


(53) 


since h, Ap and their derivatives, with respect to x, are equal to zero at \x\ —)■ 00 . Making use of the two-dimensional analogue of eq. 0: 
allows us to find: 

•I'm = y (^) ■ (54) 

Now, invoking the dimensionless variables and placing ourselves in the moving reference frame, one has: 


i/Ap( 


9C 


= -vpght I AP[-^]dU . 


Then, using the dimensionless equivalent of eq. ([^, one finds: 


-vpghl J 


■0 + c-s...») 


(55) 


(56) 


Once more, considering that C and its first derivative with respect to U are equal to zero at |f/| —)■ 00 , the first two terms from the right-hand 
side of this equation are reduced to: 


/ 


dU dU^ 


/ 


^ dU 


2 \dU 

I (C)^ 


= 0 


= 0 . 


Finally, using the definition given in eq. |2TJ, the viscously-dissipated power can be written in terms of the wave resistance: 


iVR = vr . 


(57) 


(58) 


As already mentioned, the product vr is the power required to displace the surface profile at constant speed, and it also corresponds to 
the power dissipated within the viscous thin film during such a task. 
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B Details on the two-dimensional case 

Fourier transform 


/ oo roo 

/ / ([/, Y) exp (-1 [QU + KY]) dYdU , 

-OO J —OO 

1 roo roo 

J f{Q,K)eKpii[QU + KY])dKdQ, (59) 


Wave resistance 

Starting from eq. l |35t , one considers the polar coordinates given by the relations Q = QCOs(ip) and K = psin((^), and multiplies 

both the numerator and denominator by i/V, in order to retrieve the following expression: 

roo /•27r ^3 ^ ^ (g^(f tt) 


o 2 

^ext 


dTT^V Jo 


roo r2 

Jo Jo 


cos (p2 + 1 ) 


dcpdg . 


(60) 


Additionally, if the Fourier t ran sform of the pressure field ^ under consideration is axisymmetric, a first integration of the previous equation 
over (f leads directly to eq. ( |36| >. In brief, for an axisymmetric pressure field, one should use the following formula: 


cos{^)d^p cos(^)d^ 


cos{ip)d(fi ^ 2 /* 
Jo cos{ip) + D Jo 


cos{lp) + D Jo cos{^p) + D 


= 2 


2D 




tanh ^ 


{D - l)tan 


\/l - D2 


which, considering that: 


can be re-written as: 


lim [tanh ^ (x)] = — 


f 


cos{^)d^p 


R sLtV 


cos{^) + D 

30V 

[$ (i?)] 


= 2n{l + 


iD 




L order to perform the aforementioned integration step. Moreover, a series expansion of the integrand in eq. ( |36| > near V = 0 leads to: 

i2 gdg 


which can be re-written in the form of eq. (|37| >. 

For a Dirac pressure distribution, eq. (|36|> becomes: 


Rs = 


^ r 

2ttV Jo 


A series expansion near V = 0 of the integrand yields: 


Rs = 


SL.V 


which after integration is equal to: 


27r 


Rs = 


1 - 


/■ 


2(p2 + l)' 


e +1) 


r+0{V^) } , 


(61) 


,/ v 2 + ^ 2(^2 + 1)2 


g^dg . 


gdg 


2(^?2 + l )2 


-O(V^) 




(62) 


(63) 


(64) 


Now, using the change of variables w = g + l) /V, and a series expansion of the integrand at V ^ oo, the integral in eq. ( |62| > can 
be approximated by: 


f 


1 - 


g {g‘^ + l) 


A"++ 1 )' 




^4/3 ro 

I” Jo 


f 1 _ 




-O 




4V7f 


+ 0 


which allows us to obtain eq. poll. 
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